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We discuss level schemes of small quantum-dot turnstiles and their applicability in the production of entan- 
glement in two-photon emission. Due to the large energy splitting of the single-electron levels, only one single 
electron level and one single hole level can be made resonant with the levels in the conduction band and valence 
band. This results in a model with nine distinct levels, which are split by the Coulomb interactions. We show 
that the optical selection rules are different for flat and tall cylindrically symmetric dots, and how this affects the 
quality of the entanglement generated in the decay of the biexciton state. The effect of charge carrier tunneling 
and of a resonant cavity is included in the model. 
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I. INTRODUCTION 

The constant progress in the fabrication of nanostructures 
has led to novel semiconductor devices, like quantum wires 
and quantum dots, that allow the confinement and control of 
single electrons. In quantum dots it is possible to experimen- 
tally control the tunneling of single electrons and holes. These 
systems exhibit quantum-correlations in the emission statis- 
tics [0, |[, and are very promising for future applications in 
quantum communication. 

A quantum dot that emits single photons controlled by the 
switching of a voltage is called a single-photon turnstile 
In such a system, the quantum dot is allowed to contain at 
most one single electron-hole pair, so that one photon is cre- 
ated at a time. In order to realize this, one makes use of the 
Coulomb blockade effect to suppress tunneling of a second 
electron or hole onto the dot. This implies that the system 
must be cooled to temperatures with k^T smaller than the 
Coulomb splittings. In a two-photon turnstile, two electron- 
hole pairs are created, before two successive photons are emit- 
ted. Recently, a two-photon turnstile has been proposed [Qj 
as a device to generate entangled photon pairs, which makes 
these systems very interesting. Because the Pauli principle 
allows occupation of an electronic level by at most two elec- 
trons, a two-photon turnstile can be realized without Coulomb 
blockade effects, provided that the thermal energy is smaller 
than the splitting of the single-particle levels so that tunnel- 
ing of more than two electron-hole pairs is avoided. In a 
small quantum dot, this splitting can be much larger then the 
Coulomb splittings, so that a two photon turnstile does not re- 
quire cooling in the milliKelvin regime for proper operation. 

In order to explore various possibilities for generation of 
entangled photons by two-photon turnstiles, we consider in 
this paper simple level schemes that can occur when a quan- 
tum dot is smaller than the bulk exciton size, and study which 
situations are favorable for the generation of entangled pho- 
ton pairs. In this regime, the electron and hole wave functions 
are strongly localized so that the energy separation of individ- 
ual levels is larger than the characteristic Coulomb interaction 



energies. The central idea is that, due to resonant tunneling, 
only one twofold degenerate electron level and one twofold 
degenerate hole level of the quantum dot play an active role. 
Other electron and hole levels are too remote in energy and 
may therefore be discarded in a first approximation. Within 
the resulting finite scheme of dot states we study the entan- 
glement of cascade photons from the biexciton decay, notably 
its dependence on the competition between charge carrier tun- 
neling and radiative electron-hole recombination. The effect 
of a resonant cavity will also be calculated. 



II. MODELS FOR SMALL QUANTUM DOTS 

A diagram of the semiconductor structure that we have in 
mind is shown in Figure |IJ The quantum dot is located be- 
tween P doped and N doped material. By means of a bias 
voltage V over the junction and a gate voltage <t> of an elec- 
trode near the dot, the electron level of the dot is made res- 
onant with the bottom of the conduction band of the N type 
material and the hole level of the dot with the top of the va- 
lence band of the P type material. The quasi-particle energies 
E e and Ef, are well defined as the energy of the dot with one 
excess electron, respectively hole, with respect to the neutral 
state. The bias voltage V separates the energy levels between 
the N and P sides by eV and the gate voltage shifts the elec- 
tron and hole levels by -<?<t> and e<P. The resonance condition 
for a cental dot is then 



E. - e4> = eV/2, E h + e4> = eV 12. 



(1) 
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It is energetically favorable that electrons tunnel into the dot 
when eVI2 > £ e -e4> and out of the dot when eV/2 < E e -e<P. 
The resulting level scheme is shown in Figure || and has only 
sixteen basis states, part of which are charged due to the pres- 
ence of one or two excess electrons or holes. In this scheme, 
the state with highest energy is the bi-exciton state, with two 
electrons in the upper level and two holes in the lower level. 
The opticalproperties are determined by exciton and biexci- 
ton states [pi pp. It is in the cascade decay from the biexciton 
to the ground state via a state of the one-exciton multiplet that 
an entangled photon pair may be generated. The splitting of 
the one-exciton multiplet is an effect of the Coulomb interac- 
tion between the particles [^J, which will be discussed in the 
following. 
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When one switches the gate voltage <t> or the bias voltage 
V, during a short time interval, first to a higher value and im- 
mediately thereafter to a lower value, then one promotes the 
tunneling of electrons from the N type material into the upper 
level of the dot, immediately followed by tunneling of holes 
from the P type material [gj. The system state then follows 
the path indicated by the diagonal arrows in Fig. ^. The bi- 
exciton state is produced without intermediate formation of a 
one-exciton state. Ideally the system will therefore emit a cas- 
cade of two photons, one on transition 1 and one on transition 
2, with frequencies u>\ and loi. It is important to switch the 
gate voltage <t> (and bias voltage) back to the resonant values 
(jl|) immediately after the preparation of the biexciton, in order 
to reduce the probability of electron or hole tunneling before 
the second photon is emitted. If tunneling nevertheless occurs 
before the second photon is emitted, then the transitions 3 or 
4 between the charged states may occur. Their frequencies 
CJ3 and 0J4 are in general different from u>\ and 0J2, as may 
be checked from Eq. (Q) for the energy levels. Therefore these 
photons may in principle be filtered out. In the next section we 
shall solve a master equation which includes the competition 
between tunneling and recombination processes. For small 
dots, states with more electrons or holes have non-resonant 
energies. Provided the bias voltage is not too large, the biex- 
citon state is the highest excited state and no multi-exciton 
states are formed If9fl. 
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FIG. 1 : Energy-band structure of the PIN junction, for a cross section 
through the quantum dot along the z axis. The quantum dot is a 
small cylindrical structure located in the I layer between N and P 
semiconductors. With the gate potential 4>, the electron and hole 
energies can be shifted. A bias voltage V over the junction allows 
electrons (black) and holes (white) to tunnel across the barriers with 
rate 7 and r is the photon emission rate, as indicated. 



A. Level Splitting and Photon Energies 

We consider the Coulomb interaction between electrons 
and holes as a perturbation on interaction-free levels. This 
approximation is well known in the context of quantum dots 
without holes in the electronic distribution, which are in 
(near) equilibrium [ |lO| , O], For semiconductor quantum dots 



with holes, however, localized exciton states are formed from 
electron-hole pairs. Only if the dot is smaller than the ex- 
citon size the single particle levels are well defined and the 
Coulomb interaction can be treated as a perturbation [11 12]]. 
This is the condition that we suppose to be fulfilled in the fol- 
lowing. Since we consider systems in absence of magnetic 
fields, the single-electron levels are twofold degenerate due to 
time-reversal symmetry [Q. Let the two degenerate single- 
particle states for the electron and hole level be described by 

\e)= Jdr\f){\\ Wx(7) + \l)tf 2 {7) 

\e) = fdr\f){\\ )^ 2 {r)-\i )^(r)), 

\h) = jdr\r){\] >Xitf> + li>Jfcff>). 

\h} = Jdr\7){\\) X 2{7)-\i)x\{7)), (2) 

in terms of the wave functions ipj(r) and XjO 1 ) f° r the spinor 
components. Due to spin-orbit coupling these single particle 
states are not spin eigenstates in general. The dot states that 
form the basis of the configuration space are then described 
by the occupation of the four basis states (fy. The number of 
electrons within this space thus ranges from zero up to four. 
The ground state has the hole level occupied with electrons 
and therefore is an effective quasi-particle vacuum, denoted 
with |0). Excited states of the dot are formed by means of 
creation of electrons in the higher level, and/or by creation of 
holes, i.e. by removing electrons, from the lower level. This 
gives, in addition to |0) and |e), \e), \h), \h), the further dot 
states 

\ee), \hh), \eh) \eh), \eh), \eh), 
\eeh), \eeh), \ehh), \ehh), \eehh). 

The exciton (one electron plus one hole) states may be split 
up in energy by the effective interaction between electrons and 
holes, in the form of second quantization: 

where the labels a, 0, 7, and 6 stand for e, e, h, or h. In first ap- 
proximation the antisymmetrized matrix elements V a pjs are 
those of the (screened) Coulomb interaction, but a more de- 
tailed calculation should include many-body effects. How- 
ever, whatever effects are included, in absence of external 
magnetic fields the time-reversal symmetry is always con- 
served. The one-exciton states are therefore either even (+) 
or odd (-) under time reversal and therefore given by 



\eh+) = (\eh) + \eh))/V2, \eh-) = (\eh) - \eh))/y/2, 
\eh+) = (\eh) - \eh))/V2, \eh~) = (\eh) + \eh))/V2. 



(3) 



For the same reason, the states with charge ± 1 are all twofold 
degenerate. The time-reversal symmetry implies relations like 
Veheh = Vfhfh and v eheh = v fheir The energies relative to that of 



3 



the ground state for the scheme in Fig. g may than be written 
in a closed expression as 



E = E e n e + E h n h + \{n e - l)n e V e g e e + \(n h - l)n h Vi 
+ (s- \n e n h )(V eheh + V e j w f,) + \s 2 (V eheh - V e - hen ) 



hhhh 



+ U\+s)tV eh g. + U\-s)tV. 



eheh ' 



(4) 



In this expression the symbols n e , ni, represent the number of 
electrons in the upper level and the number of holes in the 
lower level (« e ,n/, = 0, 1,2). For n e ,iih = 1, 1 we introduced 
in Eq. (Q) the notation t = ±1 for time-even or odd states and 
the symbol s has the value 5=1 for \eht) and s = -1 for \eht) 
states. For n e ,ni, =/ 1, 1, we substitute t = s = in Eq. (H). 
This splitting of the one-exciton multiplet, shown in Fig. g| 
is consistent with the phenomenological Hamiltonian used in 
refs. §0,|l5|]. 
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FIG. 2: States in a quantum dot model with one (twofold degen- 
erate) electron level and one (twofold degenerate) hole level. The 
presence of an electron is indicated by a black dot above the line; a 
hole by an open circle below the horizontal line. Besides the neutral 
states in the dot, there are states with charge ±e, ±2e where one or 
two excess electrons or holes are present. States with charge ±e are 
twofold degenerate due to time-reversal symmetry. The one-electron 
plus one-hole states are classified according to their even (+) or odd 
(-) behavior under time-reversal. Downward arrows indicate photon 
emission; diagonal arrows indicate a possible path to reach the upper 
state by subsequent tunneling of two electrons and two holes from 
the conduction and valence band into the dot. 

If the basis functions (Q) have no other quantum numbers, 
corresponding to symmetries of the system, then in principle 
all the optical transitions indicated in Fig. ^ will be present. 
This means that there is then no dark exciton state. How- 



ever, not all transition amplitudes in the decay of the biexciton 
states will be equally large. If symmetries that we consider in 
the following are fulfilled, some states will be dark. It should 
be mentioned that experimental information can be obtained 
by means of polarization measurements in combination with 
an external magnetic field [Q]. According to Eq. (Q) the sum 
of the photon energies of the transitions between the charged 
states is equal to the sum of the energies of the cascade pho- 
tons in the decay of the (neutral) biexciton state: 



(5) 



Because holes are heavier than electrons, their wave function 
are more confined leading to a stronger repulsion between two 



holes than that between two electrons, V, 



hhhh 



From this follow with Eq. (|4j) the inequalities 

U>l—U>2 > UJi - U>4 > 0. 



> Veeeg > 0. 



(6) 



The relations (Q) and (§) are expected to hold under the gen- 
eral condition that the dot is small compared to the bulk exci- 
ton size, irrespective of the shape of the dot or crystal structure 
and may therefore be helpful to analyze the emission spec- 
trum, when other information is lacking. Typical patterns of 
the emission spectra should then look similar to those plotted 
in Fig. U 

In the following paragraph we discuss the qualitatively dif- 
ferent schemes that can occur for cylindrically symmetric 
dots, and consider the implications for generation of entan- 
gled photon pairs. Total absence of spatial symmetry and neg- 
ligible spin-orbit interaction give rise to an unfavourable situ- 
ation. For then the twofold degenerate electron and hole sates 
(|^) may be written in a simpler form, without ip2 and \2 and 
with real spatial functions ip\ and \\- in that case antisym- 
metrization implies that two electrons or holes in the same 
level form a spin singlet state S = 0. The one-electron plus 
one-hole states are then a spin singlet S = 0, corresponding to 
\eh+) of Eq. (|) and a spin triplet 5=1, corresponding to \eh), 
\eh) and \eh-). Such a situation is depicted in Fig. |]d). Since 
the electric dipole operator does not act on the spin degrees of 
freedom, the S = 1 triplet will be dark and the de-exitation of 
the biexciton state proceeds only via the S = exciton state. If 
no other (spatial) symmetries are present, this is only one sin- 
gle state and consequently entanglement of the cascade pho- 
tons can never occur as that requires two different, energeti- 
cally indistinguishable paths of the biexciton decay. We there- 
fore conclude that for a possible entanglement of the cascade 
photons at least spin-orbit interaction, i.e. nonzero ip2 and \2 
in Eqs. (Q), or some spatial symmetry of the dot potential is 
required. 



B. Cylindrically Symmetric Dots Firstly the two photons must be identified as coming from 

transition 1 and 2 and must belong to the same cascade. Inco- 
herent tunneling effects between the exciton states must occur 
In order to create entanglement in photon pairs in the cas- at a slow rate 7 with respect to the lifetime of the exciton level, 
cade of Fig. ||, a number of conditions must be satisfied. Secondly, two paths in the cascade must be indistinguishable, 
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FIG. 3: Typical patterns of the photon emission spectrum for a 
(lens-shaped) quantum dot, corresponding to level scheme Qa) in the 
regime of str ong tu nneling. The plots are based the populations given 
in paragraph III A and on the relations (^|), (^j) for the photon ener- 
gies. The relative distances between the latter are arbitrarily chosen. 
The left plots shows the effect of increasing the bias potential above 
the resonance condition eV = E e + Ei, by 6, 4 and 2 times k%T. The 
right plots show, for eV = E e + E h + 2k B T, the effect of changing the 
gate potential e<I> by an amount -k%T (curve peaked at W3) to +k%T 
(curve peaked at W4) and zero (peaked at wi, as in the left figure). 



which implies that the splittings A between the intermediate 
exciton substates should be not greater than the line widths. 
If no temporal or spatial separation of the photons is possible, 
one relies on spectral separation, which leads to the conditions 
uj\ - L02 » r, r » A, and F » 7. In the general case with 
only time-reversal degeneracy of the levels, shown in Fig. [lj 
the level splittings A and lo\ - u>2 will generally be of compa- 
rable size, since they are caused by the same effective inter- 
action between the charge carriers. Identification of entangled 
pairs will also be complicated by the fact that there are four 
different routes from the biexciton to the ground state. For 
this reason we consider the case of axially symmetric quan- 
tum dots, which are often realized in experiments [^| Rh. As 
we shall see, this leads to dipole forbidden transitions, i.e. dark 
states, and degeneracy of two exciton states, so that the above 
conditions for entanglement can be satisfied. 

In an axially symmetric quantum dot, the electron and 
hole states are characterized by well defined magnetic quan- 
tum numbers ±m e and ±m/, respectively. In the most com- 
monly used semiconductor materials the conduction band cor- 
responds to si, while the valence band is a p| hole band. 
Hence the first unoccupied level in the dot is an state, while 
the highest occupied level will be a p| state. For si electrons 
and pi holes, the single-particle states (pb are of the form 



l e > = E E dr\r)(r-b\st,rnyip m (b), 

b m=-l_ J 

3 

l /z > = L E jdr\r){r-b\p\,m)Xm{b). (7) 



The summation is over the lattice sites b. The state at each 
lattice site is determined by the slowly varying amplitudes 
if> m (b), or Xm(b), and by the localized orbitals |si, m), \pl, m). 



The latter are the Wannier functions [ |l6| | which, in the tight- 
binding approximation, may be replaced by the orbitals for an 
isolated atom. The spin-components if> m (r), Xm(r) are called 
the envelope wave functions [O]. Because these are slowly 
varying with respect to the lattice, one may replace the argu- 
ment b with r in the above expressions. One finds that each 
component of the the envelope wave function is multiplied 
with a lattice periodic function, which are the Bloch states at 
the symmetry point p = 0. Projection of the states in (0) onto 
the spinor basis gives the components ip\, ip2, and \i> Xi °f 
the general expression (Q). The electron envelope wave func- 
tions ip m (r) and the hole envelope wave functions Xmff) are 
determined by solving an effective Schrodinger equation with 
an added potential U(r) that describes the position-dependent 
band edge. This leads to the confined dot states. The kinetic 
energy operator in the equation for the holes is given by the 
Luttinger Hamiltonian [jn|, ^ 

p 2 1 3 

H= 7 -^~ ? £ [73 - (73 - \' M;, 1 7", ;./,;• 

2m 9 Jim ^-J 



32 , ,,,, 9 



T i} = 3p<Pj - S u p , Jij = -(jjj +jjd - 



(8) 



The Luttinger constants 71 , 72 and 73 are dimensionless model 
parameters. The momentum operator p may be interpreted as 
a quantization of the Bloch momentum, because a plane wave 
envelope function corresponds to a Bloch wave. It is impor- 
tant to realize, however, that the physical electron position and 
momentum operators act on both the envelope wave functions 
and the orbitals in the states ([7]). 

In the envelope description of localized states in the quan- 
tum dot, one may define a total (envelope plus orbital) angular 
momentum operator as 

f=l+j = rxp+j. 

Only in the case 72 = 73, the three components of / commute 
with H [|ll , 1 7 1 . Because for InAs the two constants are nearly 
equal, this so called spherical approximation is often made. In 
the spherical approximation, a cylindrical confinement poten- 
tial U gives rise a the constant of motion mf, a spherical con- 
finement potential U results in constant / and m/. In realistic 
calculations for the case of a spherical dot [flip, one finds that 
the lowest state of the exciton as well as the lowest exciton 
and biexciton states are predominantly composed of an I = 
envelope wave function. In that case, the angular momentum 
of ground states roughly equals that of the orbital functions 
/ = ./• 

If we restrict ourselves to cylindrical dots, the confinement 
potential U is axially symmetric and the single-particle states 
have good quantum numbers mj = m e and /«/ = '«/, for the 
electron and the hole. These single-particle states will be de- 
noted as 

\e) = \m e ), \e) = \-m e ), \h) = \m h ), \h) = \-m h ), 

with positive m e , m/,. From the sixteen basis states of the level 
scheme Fig. 0, the states with an even number of electrons 
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FIG. 4: Four different schemes for cylindrical dots, a) In the case \m e - nii,\= 1, the exciton multiplet consists of a bright and a dark doublet. 



b) For the case m e 



m h 



the exciton level is split into a doublet and two singlets. One singlet is a dark state. This results in six optical 



emission frequencies, c) In case m e = m h > ^, the exciton level is split into a doublet and two singlets. Only one of the exciton states (a singlet) 
is bright, d) Level scheme for systems where spin-orbit coupling can be neglected and total spin S is a good quantum number. If no other 
symmetries are present, there is only one cascade decay path from the biexciton state. Hence, schemes c) and d) do not produce entangled 
photons. 



and of holes have total magnetic quantum number M = 0. 
One also has that V ehe i vanishes. As a consequence the pairs 
of opposite M in the one-exciton multiplet, like for example 
\eh) and \eh), are degenerate. There are qualitatively differ- 
ent schemes, shown in Fig. |j. Since electric dipole transi- 
tions occur only if \m e - mi,\< 1, we distinguish two cases: 
\m e - ni]\= 1, diagram Qa), and m e = m^, diagrams |]b) and 
Qc). The case \m e - m/,|= 1 is realized in the (lens-shaped) 
In(Ga)As/(Ar)GaAs quantum dots that have been extensively 
studied in ref. [g], There it was found that for zero external 
magnetic field the exciton states with \M\= m e + mi, = 2, which 
are formed by a m = | heavy-hole state and a m = \ elec- 
tron state, are to good approximation dark and lie below the 
\M\= 1 bright exciton states. So diagram ^) represents a re- 
alistic situation and the dark excitons with M = ±{m e + ntt,) 



may even act as intermediate stages in the formation of the 
biexciton state by tunneling of electrons and holes into the 
dot. This is also a favorable situation for the creation of en- 
tangled photons, since the two exciton states with \M\= 1 are 
degenerate, due to time-reversal symmetry, and the decay of 
the biexciton via the M = +1 and via the M = -1 exciton state 
are therefore indistinguishable, as required for entanglement. 
This only holds true, of course, if there is perfect axial sym- 
metry |jl| ^] . 

We now consider the situations sketched in diagrams |]b) 
and Qc), which represent the cases m e = = \ and m e = 
ntf, > 4 respectively. In diagram |]b) there is a doublet 
of bright states, which may allow for entanglement of pho- 
tons that are polarized in the horizontal plane. Note that the 
m/, = | state is a superposition of mi = -1,0, 1,2 states in 
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a) M = b) M = 




FIG. 5: Tunneling rates, for resonant tuning, a) The case m e ^ mi,: the bright and a dark exciton doublets are reached by equal tunneling 
probabilities, b) The case m e = m h = ^: the bright triplet and dark singlet are reached with unequal rates. Unless indicated otherwise, thin 
arrows have rate 7, thick arrows have rate 27. 



the pi hole level, while the mj, = | state is a superposition of 
mi = 0, 1, 2, 3 states. In dots elongated in the z direction, here- 
after called 'tall' dots, the ntf = \ is expected to ly below the 
nif = I state. In lens-shaped dots, the ground state has nif = | 
instead. The relevant level in a p\ hole band may therefore 
consist of the m/, = ±j states for tall cylindrical dots. The 
scheme of Fig. ^b) would also occur if, due to strain or other 
effects, the split-off pj band provides the hole states. In dia- 
gram ^p) there are two degenerate exciton states, with M = +1 
and M = -1, which are appropriate for the production of en- 
tangled photon pairs. Then one of the M = exciton states is 
dark, the other is bright. The energy of the M = states differs 
in general from that of the M = 1 states and therefore in total 
six frequencies appear in the optical spectrum. In diagram Qc) 
the exciton states with M = ±(m e + m/,) are obviously dark 
states, but also the M = time-odd exciton state \eli-) is dark. 
This follows from the time-reversal property of the dipole op- 
erator ez of the M = to M = transition. So in this case there 
is only one bright exciton state, \eh+), and therefore this situ- 
ation does not allow production of entangled photon pairs in 
the cascade decay of the biexciton. We conclude therefore that 
diagram Fig. Qb) represents a possibly favourable case for the 
production of entangled photons while a situation as depicted 
in diagram ^c) is not suitable. 

In case of a spherical quantum dot, an exceptional situa- 
tion may occur if both the electron level and the hole level are 
states with angular momentum f e = fa = \. The spherical 
symmetry then leads to a (threefold degenerate) F = 1 triplet 
and one dark F = exciton 

\eh) = \FM = 1 1), \eh) = \FM = 1-1), 
\e~h-) = -\FM = 10), \eh+) = -\FM = 00). 

In the de-excitation cascade of the biexciton state now three 
polarizations are possible for the same photon energy, which 
yield extra options for entanglement. A spherical quantum dot 



with f e or //, larger than j results in a system that can be seen 
as a combination of several systems with m e and ntf, taking all 
the possible values. The total number of states is 4^ + ^ h+l and 
the levels have a large degree of degeneracy. 

For very small dots of size comparable with the lattice con- 
stant, the crystal symmetry will be incompatible with spheri- 
cal symmetry or cylindrical symmetry. One does not expect 
that the Luttinger Hamiltonian (||) can describe this situation. 
Such small dots will fall in the class of Fig. [| Another ex- 
treme situation may arize for bound states on a single impu- 
rity atom in a further homogeneous crystal. Such a system 
resembles an ionic atom and has spherical symmetry. Our 
model applies only to the simplest situation: two levels with 
j e = j h ss i in scheme Fig. Qb), for example with a s i electron 
level and a pi hole level. 



III. STATISTICS AND ENTANGLEMENT 

A. Emission in the Strong Tunneling Limit 

An important factor that determines to what extent entan- 
gled photons will be emitted is the ratio of the spontaneous 
emission rates T in Fig. ^ and the tunneling rates 7 of the 
charge carriers. We now show how F may be experimentally 
determined in a situation of fast tunneling. For preparation 
of the biexciton the bias voltage is increased to a value where 
the tunneling rate 7 is much greater than the photon emission 
rate T, so that the electron and hole tunneling is fast compared 
to spontaneous emission. We neglect the nonradiative recom- 
bination In this regime, only thermal fluctuations can 
de-excite the system 1 10 1. When thermal energy exceeds the 
Coulomb shifts, k B T » V a /3 7< 5, the populations of the sin- 



gle particle states are independent and equal the Fermi-Dirac 
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distribution in the continuum bands: 



1 



1 



E. - e4> - eV/2 



Ph = 



1 + exp : 



E h + e<P- eV/2 



k K T 



1 + exp : 



k R T 



(9) 

The decaying levels: the biexciton, the bright exciton, and the 
two charged excitons, then have respective populations p 2 e p\, 
2p e {\ - p e )ph(l - Ph), and 2p e (l - p e )p 2 h , 2p 2 e p h (\ - p h ) for a 
flat dot (the system with \m e - mi,\= 1). For a tall dot, or an- 
other realization of the m e = mi, = j scheme, the population 
of the bright exciton is 3p e (l - p e )ph(^ - Ph) instead. Multi- 
plication of these populations with the decay rate for each of 
the levels as indicated in Fig. ||a) and b), gives the strength 
of the emission peaks. Examples of emission spectra in ther- 
mal equilibrium for the case of strong tunneling are shown 
in Fig. ||. The average emission time of a photon as a func- 
tion of temperature equals t = \l2Yp e pi l for a flat dot and is 
t = l/2(ri + T2)PePh for a tall dot. By measuring this aver- 
age one can experimentally determine Y, respectively Y\ + Ti- 
Clearly, for the preparation of a pure biexciton state, p e and 
Ph must be close to one. 



B. Correlated Photon Pairs 



The probabilities that a photon emission on transition 1 is fol- 
lowed by each of the other transitions are 



6 7 2 



Pu 2F 2 + 15r 7 + 247 2 ' * ' J 
P n =l-Pu-2P 13 , Pu = Pn- 



3r 7 + 6 7 2 



2r 2 + i5r 7 + 24 7 2 ' 



(ii) 



For tall dots, the scheme of Fig. [|b) and |]b), the emission 
probabilities on the four transitions are different, because of 
the different exciton structure. The probabilities of emission 
starting from the biexciton or starting from the exciton after 
transition 1 are in this case given by 



Pz=Pn 
P 3 = 



4 7 2 



r 2 + 9r t7 + 16 7 2 
r t7 + 4 7 2 



(12) 



r 2 + 9r t7 + 16 7 2 ' 



Pn = 



2r t7 + 4 7 2 
r 2 + 9r t7 + 16 7 2 ' 



where T t = Tj + Ti. The dependence of the probabilities in 
(flO|), ( pT| ) and ([[2]) on the ratio of the residual tunneling rate 
7 and the photon emission rate T (see Figs. Qand |]) is plotted 
in Fig. ^. It appears that the photon correlation Pu is greater 
than 90% when T > 10 7 but falls to 25% when the tunneling 
rate is much faster than the photon decay. Even if the first two 
photons are on the cascade transition 1 followed by 2, this 
does not yet guarantee entanglement. We calculate the degree 
of entanglement in the next paragraph. 



We now suppose that the quantum dot has been prepared 
in the biexciton state, so that the two-photon cascade can be 
detected. Any residual tunneling of electrons and holes can 
result in a tunneling out of the intermediate one-exciton state 
and lead to emission of a photon from another transition. At 
the resonance ([Hi, the states in the continuum levels are half 
filled, i.e. Eq. (g) gives p e = p/ t = j. Then, the unconditional 
tunneling probability of an electron into or out of the dot is 
the same and the perturbations are minimal. To obtain an an- 
alytical estimate for the relative photon emission probabilities 
and their correlations, we assume that the electron and hole 
tunneling have roughly the same rate 7. For the evaluation 
of the jump statistics, the system can be described by a clas- 
sical master equation, since only incoherent transitions occur 
[ p^ . p3| [. Since the tunneling rates do not depend on whether 
an exciton is bright or dark, the populations of the members in 
the multiplets can simply be added so that the number of rate 
equations is reduced. The net tunneling between the levels 
is indicated in Fig. ||, corresponding to the two schemes that 
can give entanglement in Fig. [|. We calculate the emission 
probabilities after preparation of the biexciton for each of the 
four transitions. The probability of a transition between a pair 
of levels is a matrix element of the inverse of the transition 
matrix, neglecting the gain terms of photon emission. For flat 
dots, the scheme of Fig. 0a) and Bp), we find the expressions 



p 2 = tL 

2r 2 + i5r 7 + 24 7 2 

P l = l-P 2 -2P 3 , P 4 = P 3 



Pi = 



2Tj + 67 2 



2r 2 + 15F7 + 247 2 ' 




FIG. 6: Photon emission probabilities for transitions 1, 2 and 3, as 
function of the ratio of the carrier tunneling rate 7 and the photon 
emission rate T of the lower transition as given by Eqs. jlo|)-([l2|). 
Thick and thin lines correspond to flat and tall dots. For tall dots 
we adopted Ti = T2 = T. Left plots: emission of the first photon 
after preparation of the biexciton; right plots: emission of the second 
photon, when the first was emitted on transition 1. The dotted line 
is the probability that no tunneling event occurs between transition 1 
and transition 2. 

The resonant tunneling rate 7, for electrons and holes, may 
be experimentally determined from the average time between 
two subsequent photon emissions in the steady state regime. 
This average time difference is for flat and for tall dots respec- 
tively given by the following two expressions: 



(10) 



_ 1 2 _ 

7 r 



2 2 
+ — + - 



sr, + 247 



97 r, 9 3r 2 + 28r t7 + 487 2 ' 



(13) 
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C. Entangled Photon Pairs 

During the switching interval that allows tunneling of two 
electrons and two holes, the system is prepared in the biex- 
citon state. This is followed by spontaneous emission of one 
photon and the system makes the transition 

\eehh) -> (\eh+)\x) + \eh-)\y))/V2, (flat dot), 
-H> (\eh+)\x) + \eh-)\y))/V2, (tall dot). 

Here \x) and |y) are orthogonal linear polarization vectors of 
the radiation field. We consider here only the case of observa- 
tion of photons emitted along the z direction. As a result, an 
entangled state between the dot and the electromagnetic field 
is formed. In order to have a degenerate doublet of bright 
one-exciton states in the the one-exciton multiplet, one needs 
an axially symmetric dot. Any asymmetry gives rise to a split- 
ting of the bright doublet of magnitude A = 2V eh f h (flat dot) or 
A = 2V e j t g h (tall dot). We consider first the case that the tunnel- 
ing rate 7 is small compared to the energy splitting A, so that 
we can neglect tunneling effects. When the system resides in 
the one-exciton state for a time f, the state will evolved into 
the state 

(\eh+)\x)e- iA ' + \eh-)\y)/V2, or {\eh+)\x) e - iA ' + \eh-)\y)/V2, 

for the respective cases of flat and tall dots. The probability 
for a waiting time t between the two photon emissions in the 
cascade equals Fe~ n , where T = Ti for tall dots. Therefore, 
the two photon density operator is given by the average 



dt Te- n (\xx)e- iA ' + \yy)) ((xx\e iA ' + (yy\) (14) 



1 



2 + 2/A/r 



\xx)(yy\+c.c. + \\xx){xx\+\\yy)(jy\. 



This expression shows that dephasing destroys the off- 
diagonal matrix element, and thereby the entanglement. The 
correlation between the polarizations \x) and |y) remain per- 
fect. We now include the residual tunneling from the one- 
exciton doublet to the four charged dot states, which occurs at 
a total rate of 47. When the system (eventually) returns to the 
bright exciton doublet, the entanglement between the dot and 
the field is destroyed and also the correlation has disappeared. 
The probability that the second photon follows the first with- 
out a tunneling event is P\ 2 = F/(F + 4j), while the second 
photon is independent of the first with probability 1 - P\ 2 - If 
we include this effect, the field density operator has a fraction 
1 - P* n that is fully mixed state, and a fraction P* 2 that is the 
average with waiting times t with probabilities re~ (r+47) '. This 
gives the two-photon density operator 



i P ( ! 

2 Vl + /A/(r + 4 7 ) 



\xx)(yy\+c.c. + \xx)(xx\+\yy)(yy\j 



+ 1(1 -P){M(x\+\y)<y\)(\x)(x\+\y){y\) 



and P = P\ 2 . It is clear from this expression that P is a mea- 
sure of the polarization correlation. The evaluation of the en- 
tanglement entropy E of such a mixed state was described in 



ref. [J24J]. After a short calculation we obtain E in terms of the 
concurrence C: 



E = -x\og 2 x-{\ 
P 

C = 



■x)log 2 (l - 
1-P 



X), 



c 2 , 



|l + iA/(r + 4 7 )| 



(15) 



When expression ( |l5| ) becomes negative, C and E are defined 
to be zero. For pure states the concurrence C gives the visibil- 
ity in two-photon interferometry [|5|] . One finds that for any 
value of A, entanglement is totally destroyed when 27 > F. 
This can be seen in the left graph of Fig. 0. As expected, the 
improvement becomes considerable when the tunneling rate 7 
is comparable with the photon emission rate F. 

One may improve on the efficiency of entangled pairs by 
detecting the photons on transitions 3 and 4. Then one can 
ignore events where the tunneling from the one-exciton level 
leads to a photon on transition 3 or 4 and only count pairs of 
photons on the cascade 1 to 2. The entanglement entropy is 
still given by Eq. (|l5|), but with P the conditional probability 
of an immediate pair 1 , 2 under the assumption that the second 
photon is of transition 2. Hence, this is P = P' 12 /Pn for flat dots 
and, provided one detects the horizontally polarized photons, 
one must put P = 3P 1 * 2 /(2Pi 2 + P{ 2 ) for tall dots. Here P 12 
is given in Eq. (|TT]) and Eq. (|l2|) the respective two types of 
level schemes. The result is plotted in the right graph of Fig. 
0. As expected, the improvement becomes considerable when 
the tunneling rate 7 is comparable with the photon emission 
rate F. 




y/r- 



y/r- 



FIG. 7: Dependence of the entanglement entropy (|l5j) on the ratio 
7/r, where for tall dots r = T[ . The upper, middle and lower curves 
are for an exciton energy splitting of A = 0, A = .2T, and A = .4r 
respectively. The left plot applies when photons from the decay of 
the charged exciton states (3 and 4 in Fig. |^) are not detected. The 
right plot applies when these photons are eliminated. Thick lines 
refer to flat dots, scheme of Fig. ^a), thin lines to tall dots, scheme of 
Fig. ^b), with axial symmetry. 



D. Quantum Dot in an Optical Microcavity 

Application of an optical microcavity (resonant with the 
lower transition 2 of the cascade), such as dielectric Bragg 
mirrors or a photonic crystal, increases the decay rate F and 
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therefore is another means to enhance the entanglement en- 
tropy. The cavity may also enhance the relative emission in a 
specific spatial direction. If, however, the cavity does not have 
two degenerate polarization modes in the xy plane, the level 
scheme of a flat dot Fig. [£|a) is perturbed and the entanglement 
is corrupted. For example, let us assume that the symmetry 
axis of the cavity is misalinged (with respect to the z axis of 
the dot) in the direction z cos 9 + (x cos <f> + y sin 0) sin 9. This 
implies that the dipole transitions corresponding to the polar- 
izations \u) = \x) cos 4> + \y) sin <j> and | v) = \y) cos <f> - \x) sin <fi 
have modified coupling constants so that the decay rates in this 
basis are T cos 2 9 and F. With quantum trajectory techniques 
fl23| , p6[ ] we obtain for the density operator 

p = fdt VGe- iH '^" 2 ( re ill '- G ' n ^'VG + (1 -P)l; 
Jo 

a=\i]xx) + \yy)) ({xx\+(yy\), 
l = \{\x){x\+\y){y\){\x){x\+\y)^\), 

which is a straitforward generalization of expression (|l4|). The 
Hamiltonian H and the decay operator G act on the state of 
photon 2 only: 

H = \xXx\A, G = | M )<M|rcos 2 9 + |v><v|r, 

and P is determined from normalization. Plots of the entan- 
glement entropy for various misalignment angles are shown in 
Fig. H Due to the complicated resonance structure of a cavity, 
the transitions 1, 3 and 4 will generally have different decay 
rates. These are preferrably smaller than the modified sponta- 
neous emission rate F of transition 2, so that Pn, P13 and P14 
are small and Pyi is nearly unity. 




in 71 

e -> 



FIG. 8: Dependence of the entanglement entropy on the relative ori- 
entation, given by the angles 6 and <f), of the cavity with respect to the 
dot. The upper, middle and lower curves are for A = .IT, A = ,2T, 
A = .4r respectively. The plots are for a constant <j> = tv/4 with a 
tunneling rate chosen at 7 = .OIT. 



the single-particle states are well separated and the Coulomb 
interaction can be treated perturbatively. We showed how 
this results in a closed level scheme with sixteen base states. 
The system seems ideal for generation of entangled photons 
on the cascade from the biexciton via the excitonic multiplet 
to the ground state. The biexciton can be prepared without 
Coulomb blockade so that low temperatures are not needed. 
For a cylindrically (but not spherically) symmetric dot, differ- 
ent combinations of the magnetic electron and hole quantum 
numbers m e , mi, give rise to the four different level schemes 
depicted in Figure || Selection rules for optical transitions 
imply that only in the first two cases, with m e - nth = ±1 or 
m e = nth = I, a degenerate intermediate level occurs in the 
cascade, which is a requirement for entanglement. Quantum 
dots in (In)GaAs/(Al)GaAs with a flat cylindrical shape have 
m e = \,m h = \ electron and hole ground states, while tall dots 
that are elongated along the symmetry axis, have m e = nth = 5 
due to restricted orbital angular momentum. Therefore, both 
level schemes, Fig. ^) and |4|b), can be realized experimen- 
tally. 

The polarization correlation and entanglement of formation 
in the photon pair may be corrupted by the following two ef- 
fects; firstly here will be a minimal residual tunneling rate 4j 
into and out of the intermediate one exciton level, which can 
effectively flip the spin of the exciton. Secondly, the Coulomb 
interaction gives rise to an exchange splitting of the exci- 
ton multiplet in dots without perfect axial symmetry, which 
causes different polarization states to dephase. The residual 
tunneling rate 7 may be obtained from the emission statis- 
tics of pairs different from the cascade 1 followed by 2 as 
given by Eq. (|l3|). The polarization correlation is found to 
be as much as 75% if 7 < 0.1T and drops only about 20% 
if 7 » r. The entanglement entropy of the two photons is 
still roughly 80% for 7 < 0.1T, provided that photons emitted 
from charged states of the dot can be eliminated. Otherwise 
is is roughly halved. The entanglement is rather insensitive to 
an energy splitting A of the (bright) exciton substates, as long 
asA<0.4r. 

Application of an optical microcavity that is resonant with 
the lower transition of the cascade, leads to increased F and 
thereby enhances the entanglement of the emitted cascade 
photons. Misalignment of the cavity axis with respect to the 
symmetry axis of the dot does not substantially decrease the 
entanglement, as long as the mismatch is less than 7f/2. 



Acknowledgments 



IV. CONCLUSIONS 

We considered realizations of a two-photon turnstile based 
on small quantum dots. In the regime of tight confinement, 



This work is part of the research program of the 'Stichting 
voor Fundamenteel Onderzoek der Materie' (FOM), which is 
financially supported by the 'Nederlandse Organisatie voor 
Wetenschappelijk Onderzoek' (NWO). 



10 



[1] P. Michler, A. Imamoglu, M.D. Mason, P.J. Carson, G.F. 
Strouse, S.K. Buratto, Nature 406, 968 (2000); C. Becher, A. 
Kiraz, P. Michler, A. Imamoglu, W.V. Schoenfeld, P.M. Petroff, 
Lidong Zhang, E. Hu, Phys. Rev. B 63, 121312 (2001); A. 
Kariz, S. Faith, C. Becher, B. Gayral, W.V. Schoenfeld, P.M. 
Petroff, L. Zhang, E. Hu, A. Imamoglu, Phys. Rev. B 65, 
161303 (2002). 

[2] D.V. Regelman, U. Mizrahi, D. Gershoni, E. Ehrenfreund, 
Phys. Rev. Lett. 87, 257401 (2001). 

[3] A. Imamoglu, Y. Yamamoto, Phys. Rev. Lett. 72, 210 (1994). 

[4] O. Benson, C. Santori, M. Pelton, Y. Yamamoto, Phys. Rev. 
Lett. 84, 2513 (2000); GS. Solomon, M. Pelton, Y. Yamamoto, 
Phys. Rev. Lett. 86, 3903 (2001); C. Santori, D. Fattal, M. Pel- 
ton, G.S. Solomon, Y. Yamamoto, Phys. Rev. B 66, 045308 
(2002). 

[5] M. Bayer, T.L. Reinecke, F. Weidner, A. Larionov, A. McDon- 
ald, A. Forchel, Phys. Rev. Lett. 86, 3168 (2001); M. Bayer, 
G. Orter, O. Stern, A. Kuther, A.A. Gorbunov, A. Forchel, P. 
Hawrylak, S. Fafard, K. Hinzer, T.L. Reinecke, S.N. Walck, 
J.P Reithmaier, F. Klopf, F. Schafer, Phys. Rev. B 65, 195315 
(2002). 

[6] E. Moreau, I. Robert, L. Manin, V. Thierry-Mieg, J.G Gerard, 
I. Abram, Phys. Rev. Lett. 87, 183601 (2001). 

[7] S. Tarucha, D.G Austing, T. Honda, R.J. van der Hage, L.P 
Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996); S. Tarucha, 
D.G. Austing, Y. Tokura, W.G van der Wiel, L.P. Kouwen- 
hoven, Phys. Rev. Lett. 84, 2485 (2000). 

[8] O. Benson, Y. Yamamoto, Phys. Rev. A 59, 4756 (1999). 

[9] P. Hawrylak, Phys. Rev. B 60, 5597 (1999). 
[10] C.W.J. Beenakker, Phys. Rev. B 44, 1646 (1991). 
[11] YZ. Hu. GieBen, N. Peyghambarian, S.W. Koch, Phys. Rev. B 
53, 4814 (1996). 



[12] A. Franceschetti. A. Zunger Phys. Rev. Lett. 78, 915 (1997). 
[13] C. Kittel, Quantum Theory of Solids, (John Wiley & Sons, New 
York, 1987). 

[14] H.W. van Kesteren, E.C. Cosman, W.A.J. A. van der Poel, C.T 
Foxon, Phys. Rev. B 41, 5283 (1190). 

[15] E. Blackwood, M.J. Snelling, R.T Harley, S.R. Andrews, 
C.T.B. Foxon, Phys. Rev. B 50, 14246 (1994). 

[16] N.W. Ashcroft, N.D. Mermin, Solid State Physics (Saunders, 
Philadelphia, 1926). 

[17] PC. Sercel, K.J. Vahala, Phys. Rev. B 42, 3690 (1990). 

[18] H. Haug, S.W. Koch, Quantum Theory of the Optical and Elec- 
tronic Properties of Semiconductors, (World Scientific, Singa- 
pore, 1990). 

[19] B. Partoens, F.M. Peeters, Phys. Rev. Lett. 84, 4433 (2000); M. 

Tadic, F.M. Peeters, K.L. Janssens, Phys. Rev. B 65, 165333 

(2002); K.L. Janssens, B. Partoens, F.M. Peeters, Phys. Rev. B 

65, 233301 (2002). 
[20] V.D. Kulakovskii, G. Bacher, R. Weigand, T. Kummell, A. 

Forchel, E. Borovitskaya, K. Leonardi, D. Hommel, Phys. Rev. 

Lett. 82, 1780 (1999). 
[21] T. Takagahara, Phys. Rev. B 62, 16840 (2000). 
[22] D.F. Walls, G.J. Milburn, Quantum Optics, (Springer- Verlag, 

Berlin 1994). 

[23] H.J. Carmichael, An Open System Approach to Quantum Op- 
tics, (Springer- Verlag, Berlin 1993). 

[24] W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[25] A.F. Abouraddy, B.E.A. Saleh, A.V Sergienko, M.C. Teich, 
Phys. Rev. A 64, 050101 (2001). 

[26] P.M. Visser, G. Nienhuis, Phys. Rev. A 52,4727 (1995). 

[27] Jian-Bai Xia, Phys. Rev. B 40, 8500 (1989). 



